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ABSTRACT 

Deep HST broad-band images taken with ACS and WFPC2 of the giant 
(~ 1000 AU diameter) dark silhouette proplyd 114-426 in the Orion Nebula 
show that this system is tilted, asymmetric, warped and photoevaporated. The 
exquisite angular resolution of ACS allows us to map the distribution of dust 
grains at the northern translucent edge of the disk, dominated by the photoe- 
vaporative flow. Using the Mie theory for standard circumstellar disk grains, 
we find evidence for a spatial gradient in grain size. The typical dust radius, 
~ 0.2 — 0.7 fim (less than what reported by previous studies) becomes smaller as 
the distance from the disk center increases, consistent with the expectations for 
the dynamic of dust entrained in a gaseous photoevaporative wind. Our analysis 
of the disk morphology and location within the nebula indicates that this system 
is photoevaporated by the diffuse radiation field of the Orion Nebula, while be- 
ing shielded from the radiation coming directly from the central Trapezium stars. 
We estimate the mass-loss rate from the disk surface and the time-scale for total 
disk dissipation, which turns out to be of the order of 10 4 yr. Such a short time, 
of the order of 1/100 of the cluster age, indicates that this system is seen on the 
verge of destruction. This is compatible with the exceptional nature of the disk, 
namely its combination of huge size and low mass. Finally, we briefly discuss 
the viability of possible mechanisms that may lead to the peculiar morphology 
of this system: external UV flux, binary star and past close encounter. 
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Introduction 



Early Hubble Space Telescope (HST) images of circumstellar disks in the Orion Nebula 
Cluster have provided the first direct evidence that protoplanetary disks in rich young stellar 



clusters may be photo-evaporated by UV radiation of nearby massive stars (O'Dell & Wen 



1994). The typical signature of photo-evaporation is a bright ionized cusp pointing to the 



brightest stars of the cluster, ^ 1 Ori-C in particular. The cusps always surround a central 
low-mass star, that may occasionally be encircled by a compact dark disk seen in silhouette 
against the nebular background, or against the cusp itself. The most recent atlas of pro- 
toplanetary disks ("proplyds") in the Orion Nebula (Ricci et al. 2008) counts 178 of these 



comet-shaped sources. The same atlas also counts other 36 disks which do not show evi- 
dence of photoevaporation, being seen as pure dark silhouettes either in absorption against 
the nebular background (28 sources) or as dark lanes at the midplane of some diffuse bipolar 
emission (8 sources). Disk photo-evaporation naturally raises the issue of disk survival in 
the harsh environment of young massive clusters (see e.g. IQ'Dell et al.|2008] ). The frequency 
and spatial distribution of disks with and without evidence of photo-evaporation within the 
same cluster provides an excellent opportunity to investigate how environmental factors may 
affect the path to planet formation. 



The basic physical process at the origin of disk photo-evaporation in massive clusters 



has been clarified by Johnstone et al. (1998). The FUV (6 eV < hu < 13.6 eV; 2000A> A > 



912A) radiation reaches nearly undisturbed the disk surface, creating a photodissociation 
region (PDR) with temperature T ~ 100 — 1000 K. Depending on the radial distance from 
the ionizing stars, the heated disk surface may reach thermal velocity higher than the escape 
velocity. The resulting neutral outflow, largely composed by molecular H 2 gas dragging a 



non negligible amount of small dust grains (Owen et al. 2011) soon encounters the EUV 



(hu > 13.6 eV; A < 912A) radiation coming directly from the ionizing stars, creating a 
ionization front at T ~ 10 4 K, i.e. the bright ionized cusp. This model can explain both 
the general morphology of photoevaporated disks as well as more peculiar features, like the 
presence of [OI] emission from photodissociation of OH radicals produced by H 2 heating and 
dissociation observed in one of the most prominent objects, 182-413 p|( |Storzer fc Hollenbach 



1999 Bally, O'Dell & McCaughrean 2000). 



In this paper we shall adopt the coordinate-based naming convention of 



O'Dell et al. 



(2008) 
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It is generally assumed that disk photo-evaporation in rich clusters, as revealed by the 
presence of photo-ionized cusps, can occur only if the disk is directly exposed to the EUV 
flux from the massive stars, consistently with the fact that the majority of photo-evaporated 
disks in the Trapezium cluster lies in the immediate surroundings of 6 |1 0ri-C Dark silhouette 
disks, lacking any evidence of photo-ionization, are believed to lie within the foreground Veil 



of the Orion Nebula, shielded by EUV radiation (McCaughrean & O'Dell 1996). However, 



dark silhouette disks can also be photo-evaporated by FUV radiation (Adams et al. 2004) 



although the products of photo-dissociation in these systems have never been seen (O'Dell 



et al. 2008). This because the parameter controlling photo-evaporation is the intensity of 
the FUV flux, which can be expressed in units of Habing Flux (the estimate average flux in 
the local interstellar medium), defined as Go = 1.6 x 10~ 3 erg cm -2 s _1 . In the inner core of 
the Trapezium, the FUV flux is elevated, reaching G = 10 5 Go at the Orion I 
, about 110"or 0.22 pc south of 6 )1 Ori-C (assuming a distance d = 4 



et al. 



et al. 



2002 



(Robberto 


pc 


Menten 



2007) Sources in the outskirt of the region, shielded from the EUV photons emitted 



e.g. by 6' 1 Ori-C, may still receive scattered EUV and, especially, FUV radiation from the 
nebular environment. 

The most spectacular case of dark silhouette in Orion is the well know 114-426 disk 



(McCaughrean et al. 1998). With an apparent diameter of over 2", which makes it visible 
also from the ground, this disk spans about ~ 1000 AU, 10 times the size of our Solar System 
or of other typical photo-ionized disks in Orion. It is seen nearly edge-on, the presence of 
the central star being noticeable through scattered polar emission, and is considered to be 



in a relatively quiescent, non photo-evaporated phase (Throop et al. 2001). McCaughrean 



et al. (1998) estimate for the central star an initrinisic brighness K ~ 9.5 m , corresponding 



to a 1.5M star at the distance and age of the Orion cluster, but the uncertainites are large. 



Using HST/NICMOS high-resolution near-infrared images of 114-426, McCaughrean et al. 
( |1998[ ) find that the major axis of the disk in the Paa at 1.87 /im line is ~20 % smaller than at 
0.6 /im. From this, they deduce a minimum disk mass of 10 M e , and possibly > 5 x 10" 4 M . 



On the other hand, 114-426 has not been detected by the Submillimeter Array survey (Mann 



& Williams 2010) at either 880 /mi or 1330 /mi. This non-detection places an upper mass 
limit of 1.2 x 1O~ 2 M (0.5 - 12 M Jup ), confining the disk mass into a range of two orders of 
magnitude. 

This object has attracted interest because the variation of the disk size with wavelength 
has allowed to probe the properties of the disk grains, exploiting the wavelength depen- 
dence of the extinction of the bright nebular background through the translucent disk edges 



2 while Menten et al. (2007) provide the most precise estimate, O'Dell & Henney (2008) consider a distance 
d = 430±20 pc to be more accurate. Our estimates do not depend on the exact value of the assumed distance. 
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(Throop et al. 2001). If the grain sizes were much smaller than the observed wavelengths, 
as in the Rayleigh scattering limit, shorter wavelengths would be scattered much more ef- 
ficiently than longer ones: in this case the disk would appear smaller at long wavelengths. 
Viceversa, if the grain sizes are of the same order of the wavelength of observation, the disk 
size is expected to remain nearly constant with wavelength ("gray" scatter). Using obser- 
vations of 114-426 in the Ha (0.656 /jm; HST/WFPC2), Paa (1.87 /im; HST/NICMOS1) 
and Bra (4.05 //m; Keck/NIRSPEC), Shuping et al.| ((2003) found evidence for chromatic 



extinction at the northern translucent edge of the disk, deriving a typical grains diameter 
> 1.9/mi, but not ^> 4 /mi. This represents one of the first evidences for grain growth in 



protoplanetary disks, the initial step in the process of planet formation (Natta et al. 2007) 



In this paper we use the most recent HST/ACS images of the 114-426, together with 
a HST/WFPC2 image taken in the U-band, to clarify the disk structure and its physical 
status. Contrary to what is commonly assumed, we show that the disk appears strongly 
affected by the action of the environment and is photoevaporated. Its unusually large size 
and morphology (warped and asymmetric with respect to the central star) may result from 
photoevaporation. The translucent northern edge of the disk appears originated by the 
outflow of neutral material rather than the truncated edge of an undisturbed disk. We 
map the variation of extinction in the translucent region, showing how the typical grain size 
variation is compatible with an evaporative flow. This allows us to estimate the mass-loss 
rate from the disk surface and the disk lifetime. We also show that the disk is enshrouded 
by a large scale (~ 0.03 pc) cocoon of dust with a bright ridge facing the side opposite to 
the cluster center. This suggests that the disk photoevaporation is driven by the diffuse 
radiation field rather than by direct irradiation from the massive cluster stars. 

In Section [2] we briefly present our data; in Section [3] we illustrate the geometry of 
the source, moving from the inner structure to the most extended features and the nebular 
surroundings; in Section [4] we analyze the grain size composition and distribution in the 
photoevaporated flow; in Section [5] we derive the main physical parameters of the mass loss 
and discuss the possible heating sources. Our conclusions are presented in Section [6} 



2. Observations 

HST/ACS-WFC and HST/WFPC2 images of 114-426 have been collected as part of the 
HST/Treasury Program on the Orion Nebula (GO-10246, PI M. Robberto). In this paper 
we use the data taken in six broad-band filters: F336W with WFPC2, and F435W, F555W, 
F658N, F775W, F850LP with ACS (Fig. [T|. The passbands are approximately equivalent 
to the standard U, B, V, Ha, I and z photometric system, respectively, and for simplicity 



we will equivalently use both naming conventions. 



The instrument configuration, observation and data reduction strategy are presented 
elsewhere (Robberto et al. 2012). Here we remind that the ACS images have been drizzled 
to a scale 50 mas/pixels, corresponding to the nominal scale of the instrument, whereas the 
WFPC2 images have been coadded preserving the original 10 mas/pixel scale of the PC 
chips. The FWPC2 image has been resampled and rotated to align it to the ACS images. 
A column of bad pixels in the WFPC2 image survived image processing and is visible in the 
F336W image. However, it's location and extent does not affect our analysis. 



3. Morphology of 114-426 

3.1. Disk Structure 

In Fig. [T] we present our six images ordered in order of increasing wavelength. The 
scattered light from the central star, leaking through the polar lobes of the disk, becomes 
more prominent with wavelength. The western bright lobe is always brighter than the eastern 
one, indicating that the disk, seen nearly edge-on, is tilted with the western side facing us. 
The morphology also shows clear signs of multiple asymmetries. First, the lobes are not 
symmetric with respect to the disk's polar axis, being displaced from the center and closer 
to the southern edge. Second, the lobes are not on the same vertical axis: a line joining the 
centers of the two lobes would clearly appear tilted with respect to the disk axis, with an 
inclination of ~ 45°. 

To further investigate the structure of this system, we have subtracted from the F850LP 
image the F658N image, properly scaled, to enhance the image contrast by exploiting the 
cromatic variation of dust opacity (Fig. [2]). The disk appears warped, with an "equatorial" 
dark lane tracing in silhouette the asymmetry of the bright lobes. The southern edge of the 
disk is remarkably sharp and thick, at a projected distance of ~ 370 AU from the mid-point 
between the two lobes, which may be assumed as the position of the central sta . We will 
take Td = 400 AU as the nominal radius of the disk. 

In what concerns the northern edge of the disk, by stretching the gray levels of our 
F555W image we can resolve it in a spur-shaped structure (Fig. [3j protruding beyond the 
darkest part of the disk. If we regard the spur as an outer appendix to the main disk, 
thus limiting the disk extension at the northern side, the symmetry of the warped disk with 
respect to the mid-point between the two lobes (i.e. the position of the central star) is nicely 
recovered. We believe that the spur most probably represents an outflow of photoevaporated 



material. It corresponds to the "translucent" part of the disk where Throop et al. (2001) 
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F336W ACS F435W ACS F555W 




Fig. 1. — HST images of the protoplanetary disk 114-426 in our six photometric filters, 
equivalent to the standard U, B, V, Ha, I and z passbands. 




Fig. 2. — Difference of images of 114-426 taken in F658N and F850LP filters, showing the 
warped geometry of the system. 
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and Shuping et al. (2003) have measured the dust grain size. In Figure Else shows a cartoon 



of the system summarizing the main features: tilt, asymmetry and warp. 



3.2. Dark filaments 

Observing on a larger scale the Ha image of the source, one can detect further signspots 
of ouflow of neutral material. Fig. [5] shows both a color composite image and a contranst 
enhanced V-B image of the target, stretched to enhance the contrast against the nebular 
background. Both images reveals a remarkable dark "tail" on the eastern side, characterized 
in the innermost region by twisted filaments that can be traced up to distances of several 
disk diameters. The darker filament departs from the southern part of the disk, just to the 
south of the eastern lobe, and can be directly revealed by a contour plot of the disk intensity 
levels (Fig. [6]). The same figure also shows a minor profusion on the north-western side of 
the disk, just to the north of the bright lobe. Notice that these two dark protrusions appear 
simmetric with respect to the disk center and antisymmetric with respect to the bright lobes. 



3.3. 114-426 in its Environment 

Figure [7] shows our source with its immediate surroundings in a contrast enhanced Ha- 
V image. An extended filamentary shell emitting in Ha line encircles the western side of the 
disk. The shell radius has a projected distance of about 2.5 disk diameters, ~ 0.01 pc. The 
filaments can be traced, at much lower brightness levels, to the eastern side of the disk, up to 
a projected distance of 0.03 pc. The tenuous bright rim delineates a characteristic "foot-like" 
structure which seems to represent a region of enhanced density externally illuminated by 
the diffuse radiation field of the nebula. 

Remarkably, the two closest sources immediately to the east of 114-426, 117-421 and 



121-434 (Ricci et al. 2008), are resolved in our ACS images in two compact photo-ionized 
proplyds, pointing in the direction of ^ 1 Ori-C, lying about 100" to the north-east (~0.2 pc 
projected distance). The small size of the two dark disks, unresolved in our ACS images and 
in any case less than ~ 50 AU radius (bounded by the ionized cups), is in striking contrast 
with the huge dimensions of our target, about 20 times larger in size and without any sign 
of direct influence from 6 |1 Ori-C. 



In conclusion, there are various elements of circumstantial evidence all pointing in the 
directions of 114-426 being photoevaporated, and we uphold that this is indeed the case. 




Fig. 3. — Image of 114-426 in F555W filter. The upper edge of the disk presents a "spur- 
shaped" protrusion with a projected dimension of ~ 250 AU. 




Fig. 4. — Three graphic models of the tilted (a), warped (b) and photoevaporated shape of 
the disk. 
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Fig. 5. — Top: Color composite image of 114-426 in the F435W (blue), F555W (green) and 
F658N (Ha, red) filters; Bottom: difference of the F555W-F435W images, zoomed-in on 
the target. Both images show evidence of dark filamentary "tails" of photoeaporated dust 
escaping from the disk. The disk size, about 2"diameter, provides the images scale. 
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Fig. 6. — Countour plot of 114-426 in the F658N (Ha) filter, showing the antisymmetry of 
the dark filaments with respect to the location of the bright lobes. 
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Before envisioning a possible scenario for this unique object, it is appropriate to reconsider, 
on the basis of our data, the physical properties of the dust grains. This is the subject of 
the next section. 



4. Dust particles sizing 



As shown by Throop et al. (2001) and Shuping et al. (2003), it is possible to investigate 



the properties of dust particles of 114-426 via the wavelength dependence of the extinction of 
background light in partially-opaque areas of the disk. Following the same type of analysis, 
we are now able to exploit the spatial resolution of our ACS images to map for the first time 
the particle size in different positions of the partially-opaque part of the disk. 

We use to this purpose the classic Mie theory, a fully analytical model which derives, for 
the special case of spherical particles, an exact series expansion of the scattered field (see for 
example Van de Hulst|1957 ). Three parameters are needed for computing the cross section of 
dust: the radius a, the real part of the refractive index Re(n), and the imaginary part of the 
refractive index Im(n), the latter two quantities being dependent on the wavelength of the 
incoming light. Both the absorption (C a b s ) and the scattering (C sca ) cross sections are thus 
determined, and therefore the extinction cross section given by C ext = C a b s + C sca . Generally 
speaking, the dominant role is played by the radius a, while the refractive index components 
have a relatively minor effect. This fundamental property allows performing a measurement 
of the grain size from the shape of the cross section curve as a function of wavelength. The 
refractive index of the grains can be estimated given the dust grain composition, which is 



typically considered to be a compound of different chemical species (e.g. Pollack et al. 2004 



Semenov et al.|20"03 ). The Bruggeman mixing rule (Bruggeman 1935 ) provides a generic form 
that expresses the dielectric function of the heterogeneous particle, eo,e//, as a combination 
of the dielectric functions of the single species. For a iV-components mixture it is: 



iV 



e 0,i 



6 0,eff 



e 0) i + 2e , e // 



where 6q j and /j are the dielectric function and the volume density of the i — th species 



(X)i=i fi = 1)> respectively For the dust particles of circumstellar disk we consider porous 
composite grains made of astronomical silicates (/ = 0.07), carbonaceous materials (/ = 
0.21), water ice (/ = 0.42) and vacuum (/ = 0.30). This is an approximated version of the 
Pollack et al"] ( |2004| ) model (see |Ricci et aL]|2010[ ) . The dielectric function is derived from the 



complex refraction index, £j 



We have assumed for the Silicates the refraction index 
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13 



vs. the wavelength tabulated by Weingartner & Draine (2001) for the Silicates, by Zubko 



et al. (1996) for the amorphous Carbonacious material and by [Warren (1984) for the Water 
Ice. For vacuum, eo }Va c = 1- By solving the implicit Equation [l] for e e // we derive for the 
real part of the refractive index Re(n) ~ 1.3 and an imaginary part, Im(n) ~ 0.03. Having 
derived the effective refraction index (or dielectric function), we derive the absorption and 
scattering cross sections, and therefore the extinction, from the Mie theory. 

The shape of the dust grains deserves some discussion. Real dust grains are not spherical 
but likely present complex internal structures due to the aggregation of smaller interstellar 



grains. The most recent studies (Min et al. 2006 Borghese et al. 2007) indicate that the 



scattering properties of particles of different shape and sub-structure are not dramatically 
different from those derived using the Mie theory for spheres. Assuming that the dust 
particles are perfect spheres does not introduce uncertainties larger than a 30%, the size 
returned by the classic Mie theory being typically underestimated. 

Also the size distribution plays a role in determining the extinction. In this paper we 
make the simplifying assumption of treating the dust grains as if they were monodisperse, 
just one radius, in each pixel. In a separate paper the issue of the size distribution will be 
treated more in detail. 

To compare theoretical predictions with observations, we exploit the extinction efficiency 
factor, Q ex t, namely the cross section divided by the geometrical cross section G = ira 2 . This 
is a parameter commonly used in scattering problems, and can easily be obtained on the 
basis of the Mie expansion series as a function of the reduced size parameter x = ka, where 
k = 2tt/\. The extinction curve can then be obtained at any A for a given size a, and then 
fitted to the data. The extinction data have been evaluated as the ratio of the average flux 
(counts/s) of the nebular background, F\ )backl and the flux (counts/s) in a single pixel of 
the translucent zone^J F^ pix . In particular F\^ack has been evaluated in an elliptic ring, all 
around the disk, close enough to be rapresentative of the effective nebular background of 
114-426. More precisely we define Q ex t, given by the observational data, as 



Q 



ln( 



( .it 



(2) 



The surface density X disappears when Q ex t is normalized to some reference wavelength, 
for example the V-band (F555W). The plots in Fig. [8] show the data for three representative 
pixels and the corresponding best fit curves, providing radii equal to 0.60 /xm, 0.50 /im and 



3 Thc A index labels the broad-band filter. 
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0.36 fjjn.. The attenuation in the V-band (F555W) varies from F V fi ack /F V)Pix = 4.0 (pixel A) 
down to Fy t back I Fy, P ix — 1-17 (pixel C). We estimated the error bars accounting for both the 
Poisson statistics (on F\^ ac k and F\ tPix ) and non uniformity of the background (for F\^ ac k). 
Similar plots have been obtained for every pixel in the translucent zone, allowing us to map 
the distribution of dust grains size, also shown in Fig. [8] From the outmost pixels (blue) 
to the innermost pixels (light yellow), the radii increase from 0.2 jttm to 0.6 /im. The inner 
parts of the disk cannot be studied with this approach since the optical thickness increases 
and the scattered light from the lobes dominates. 

The absolute values of the radius a we have derived depend on the assumed grain 
composition. For a generic grain of refraction index n, the scattered field depends on the 
product xn = 2nan/X: an increase of the refraction index n compensated by a similar 
decrease of the grain radius a provides the same result. Therefore, the uncertainties on the 
assumed refraction index affect only the absolute value of the grain size, not their ratios. 
The evidence for a spatial variation of the grain sizes at the disk edge appears therefore 
strong. 

Our analysis, based on Eq. (2), does not account for a scattering term. Heuristically, 
we not expect scattering to be dominant. One the one hand, it is reasonable to assume that 
our dust particles are non-spherical, randomly oriented, heterogeneous and with variable 
refraction index. In this case, the oscillations of the individual cross sections predicted by 
the Mie theory are smeared out and no large-scale chromatic patterns should be expected. 
Also, scattering by typical disk grains is highly non-isotropic, as Mie theory predicts 90% of 
the scattered light is emitted in the forward hemisphere. Since our source is mostly back- 
illuminated, wave propagation toward us is dominated by the coherent superposition of the 
background light with the forward scattered wave, i.e. extinction. Third, even assuming 
isotropic scattering, it may be shown that the color effect induced by the scattered radiation 
produces an apparent modulation of the reddening law (W. Henney, private communication) 
that seems unable to explain the observed range of attenuation, and therefore of grain size, 
shown by our data. 



5. Discussion 
5.1. Grain size variation in the photoevaporative wind 

Some evidence for the spatial variation of the size of dust grains has been obtained by 



recent interferometric observations of young disks in the sub-millimeter (Banzatti et al. 2011 



Guilloteau et al. 2011 ). The general result found so far is that larger grains are found in the 
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Fig. 8. — Size particles distribution in the photoevaporated area at the NE edge of 114- 
426. Different colors represent distinct mean radius a. The three plots show the mea- 
sured absorption efficiency and corresponding extinction curves for the labeled pixels. 
The image is pixelated according to the ACS/WFC image scale, 1 pixel=50 mas=21 AU 
at d=414 pc. The coordinates of Pixel A are RA(2000.0)=5:35:11.345; DEC(2000.0)=- 
5:24:25.70, from the ACS images of the Orion Nebula obtained by the HST Treasury Program 
(file orionJist_acs_StripOLf658n_vl.fits on http : / / archive .stsci.edu/ prepds /prion/ \ . 
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disk regions which are closer to the central star, and this is the same trend found for the 
114-426 disk. In those cases, however, the dust probed by those observations is located inside 
the disk, and not in a photoevaporative wind as for the translucent north side of 114-426. 

The result presented in this paper is in general agreement with the expectations for 
the dynamics of dust grains entrained in a gaseous photoevaporative wind. The drag force 
experienced by a dust grain is proportional to a 2 (Takeuchi et al. 2005), and therefore 
the acceleration scales with a -1 . This means that the wind selectively entrains grains with 
different sizes, with smaller grains being more easily dragged along the wind. As shown 



by Owen et al. (2011), the drag force on small grains can dominate over the stellar gravity 
allowing them to be carried to very large distances from the central star (several hundreds 
of AU). This can explain, at least in a qualitative way, the radial gradient of the grain sizes 
obtained by our analysis over a spatial scale of 200 — 300 AU (Fig. ??). 

However, the Owen et al. models treat the case of a photoevaporative wind driven by 
the ionizing radiation from the central star, whereas in the case of 114-426 the observed 
asymmetries indicate that the origin of the photoevaporative wind is external to the system, 
i.e. due to the diffuse FUV radiation field of the nebula (Section 3). This has an obvious effect 
on the geometry predicted for the photoevaporation of the disk. A model which calculates 
the photoevaporative wind induced by a diffuse non-ionizing radiation field is needed to 
obtain a more quantitative comparison with the data. 



5.2. A scenario for 114-426 

To build a possible scenario for 114-426, we begin by observing that the disk surface 
appears shielded by the EUV flux from ^ 1 Ori-C, unlike the nearby proplyds 117-421 and 
121-434 (Fig. [7]). This can be explained in two ways. On the one hand, 114-426 may lie 
much further away from the Trapezium than the other two proplyds, possibly within the 
foreground Veil of the Orion Nebula, shielded by the EUV photons from ^ 1 Ori-C by some 
extended ionization front, possibly at the inward-facing side of the Veil. In this case, it seems 
difficult to explain the existence of the bright filamentary structure encircling our source on 
the West side (Fig. [7]), which is basically the direction opposite to 6 1 Ox\-C This filament, 
mostly visible in Ha, most likely represents a ionization front produced by diffuse EUV 
radiation coming from the outer regions of the nebula, absorbed by a density enhancement 
apparently centered on our target. At the same time, there is no evidence of a comparably 
strong ionization front on the East side, facing 6 |1 0ri-C In order to be viable, the "within 
the Veil" scenario seems to require a rather peculiar combination of geometry, veil density 
and radiation field. 
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The second possibility is that the disk lies in the outer parts of the nebula in a location 
shielded by the EUV radiation directly coming from 6 |1 0ri-C It is intriguing, in this respect, 
to notice that between 114-426 and the Trapezium lies, projected on the plane of the sky, 
the high-density region hosting Orion-S, a secondary center of star formation in the Orion 



Nebula (Fig. ph. According to most recent 3-D model of the Orion Nebula (O'Dell et al. 



2009 O'Dell & Harris 2010), the Orion-S feature southwest of the Trapezium is an optically 
thick molecular cloud within the body of ionized gas. The Orion-S Cloud that may therefore 
prevent EUV photons from the Trapezium stars to reach 114-426. The two proplyds 121-434 
and 117-421 would therefore lie in the foreground, far enough from the Orion-S Cloud to see 
directly the Trapezium stars. Alternatively, they could be located deeper inside the Nebula, 
exposed to the EUV flux passing "under" the Orion-S Cloud. 

We propose that the influence of the FUV photons mostly coming from the Western 
regions of the Orion Nebula causes 114-426 to photoevaporate. The contribution from other 
stellar sources to the East of 114-426 is neglible, in particular we estimate that the total flux 
from closest AO star KS Ori (Par 1685) at ~ 3 arcmin distance is about 1% of the diffuse 
FUV field (see the next section). The neutral outflow from the disk surface intercepts the 
weak diffuse field of EUV photons at relatively large distances from the disk; still, the ionized 
filamentary structure may well represent the analog of the bright ionized cusps that identify 
common proplyds. The pressure gradient pushes the neutral material further to the East, 
originating "foot-like" structure traced by a weak ionization front at its boundaries. Also 
the photoionized proplyds in the vicinity of theta 1 Ori-C show dusty tails downstream of 
the ionizing radiation. In summary, 114-426 can be regarded as an extreme case of proplyd 
photoevaporated by the diffuse radiation field. 



5.3. Mass loss from 114-426 



The mass loss rate, M, for large disks that evaporate in an external FUV field is given 



by (Adams et al. 2004): 



M = A-KF(n)a F }j V c s r d . 



(3) 



where J 7 is the fraction of solid angle subtended by the outflow, fi is the mean molecular 
mass of the gas, <Jfuv is the absorption cross section for dust grains to FUV radiation, c s is 



the sound speed of the flow, and is the disk radius. With Adams et al. (2004), we assume 
T ~ 1. Assuming a typical dust grain radius a = 0.5 /im, as derived in Section |1J and a 
gas-to-dust mass ratio of 100 : 1 we derive <Jfuv — 8.4 X 10~ 22 cm 2 , in line with the typical 
value for PDR regions (Storzer & Hollenbach 1999). For the outer radius we use ~ 400 



AU, concentrating our analysis on the northern translucent region. The sound speed, i.e. 
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Fig. 9. — Position 



of 114-426 with respect to the Trapezium, Orion-S and LL Ori 
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the temperature at the disk surface, is controlled by the amount of FUV illuminating the 
disk which can be estimated as follows. 

Both ioinizing (EUV) and non-ionizing (FUV) fluxes emitted by the Trapezium stars 
are scattered within the HII region. The EUV radiation, carrying about 2/3 of the original 
ionizing energy (assuming Menzell's case B recombination), is scattered by recombination 



processes, whereas FUV radiation is scattered only by dust particles (see e.g. Robberto et al. 



2002). Most of the scattering occurs at the edges of the HII region so that the diffuse field 
can be assumed to be constant within the cavity. A disk shielded from the EUV radiation 
and from direct illumination from the cluster center will only receive a flux of scattered FUV 
photons given by 

F ,. ff „ a _Q . LfUV U) 

assuming that 1/2 of the radiation is redirected inward by the walls of the cavity, a scattering 
efficiency Q sca ~ 1 appropriate for large grains, a flux Fpuv of non-ionizing photons emitted 
by the central O stars (mostly 6 |1 Ori-C), and a typical radius of the cavity Rhii = 0.3 pc, 
the distance between ^ 1 Ori-C and the ionization front ( O'Dell et al.||2008 ). We assumt^jthat 



^ 1 Ori-C emits 50% of its luminosity as FUV photons, L EUV - 1.5 x 10 49 photons s 1 and an 
effective wavelength Xpuv = 1200 A. A shielded disk will then receive a diffuse FUV flux 

? diff 



F diff ~ 2 • 10 3 G . 



Solving the equilibrium equation for the radiative heating and cooling of dust grains 
(e.g. |Draine |[20lT| Ch24), assuming; the above FUV flux and the refraction index of our het- 



erogeneous grain, we derive a dust temperature of 25 K, appropriate for a Photo Dissociation 



Region. For the temperature of the gas components, we can refer to Fig. 2(b) of Adams et al. 



(2004), which reports the temperature profiles calculated from the PDR code for an external 
FUV radiation field of 3000 Go as a function of visual extinction Ay, for different values 
of the number density. At Ay ~ 1 the gas temperature is in the range T g ~ 100 — 300 K, 
corresponding to a sound speed c s ~ 0.6 — 1.0km s _1 from the disk. 

Using Equation (3) we then derive for the mass loss rate of 114-426 

M~3-5xl0~ 7 — (5) 

yr 

where the spread is mainly dependent on the uncertainty on the sound speed. 



4 6 |1 0ri-C is historically flagged with a spectral type 06.5pe. However, recent observations indicate that 
this peculiar system is a binary with a 05.5 primary (M=34.0M Q ) and an 09.5 secondary (M=15.5 MQ 



Kraus et al.|2007(@ ). We therefore adopt the Lyman continuum luminosity of a ZAMS 05.5 star, according 



to 



Panagia (19731. Our results to not depend on this assumption. 
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5.4. Disk lifetime and unicity 



The large disk size, combined with the modest upper limit for the disk mass, make 
114-426 a rather unique object. Analyzing the thermal continuum (880 /xm) emission from 
circumstellar disks in the ~ 1 Myr old Ophiucus star-forming region, Andrews et al. (2010) 



have found a strong correlation between the disk mass and radius, larger disks being more 
massive. A comparison of our values for the disk radius and the mass upper limit with the 



valuer reported by Andrews et al. (2010) shows that our source does not follow the radius- 



mass relation traced by the other sources, being way too large for its mass. If the relation 
radius-mass is an imprint of the initial conditions inherited by the disk at the formation 



epoch, as suggested by Andrews et al. (2010), then 114-426 is in a diverse evolutionary stage 



in which this dependence seems to be lost. 

Assuming our range of estimated mass loss rates and the current upper limit for the disk 
mass, 1.2 x 1O~ 2 M , the disk's life time ranges between 2.2 x 10 4 yr to 3.7 x 10 4 yr. This is 
clearly a very short time scale, but it may be appropriate for a rather unique object like 114- 
426. At the current evaporation rate, 114-426 will be entirely dissipated in a few ~ 10 4 yr. 
It may well have formed 1 Myr ago (about the age of the ONC) with an average disk mass 
of ~ 0.1 — 0.5 M and be reduced to the current mass maintaining level of photoevaporation 
nearly constant at the measured level. Overall, the picture seems consistent with 114-426 
being a relatively low mass disk inflated and photoevaporated by the diffuse field of FUV 
radiation within the nebula, seen in the very last phases of its evolution. 



5.5. Disk warping in the ONC 

In this last section we briefly discuss the disk morphology. First, we note that the 
keplerian orbital period at the outer disk edge 

/ 47r 2 r 3 

Pout = V gm1 ^ 8000yT (6) 

corresponds, for dust grains moving at lkm/s, to a traveled distance of about 1700 AU, or 
4 disk radii. This is nicely compatible with the spatial scale of variation of the spiraling 
dust filaments on the eastern side of the disk, which may therefore be related to periodic 
enhancement of the disk mass loss rate syncronous with the passage, or appearance, of disk 
asymmetries. 

Concerning the disk asymmetrywhile the diffuse FUV radiation represents the main 
mechanism of destruction of 114-426, it cannot easily account for the extremely warped 



disk geometry unless one assumes it is highly non-isotropic. A directional radiation flux 
impacting the surface of a flared disk can generate a non-uniform increase of the surface 



temperature and modulate the flaring angle Robberto et al. (2002), or just push the disk 



mid-plane out of its equilibrium position, depending on the disk angular velocity and surface 
density (S. Lubow, private communication). Since the disk is rotating, the development of 
a persistent warp will depend on factors like the cooling time scale, the turbulence regime, 
or any possible temporal variations of the radiation intensity, all basically unknown at the 
moment. 

A second possibility is that the disk hosts a binary star, in which case the warping would 
result from a a misalignement between the disk plane and the plane of the binary orbit. 



Very little is known about the nature of the central star (see McCaughrean et al.||1998 for a 
discussion). An intriguing aspect is that the colors of the two lobes are quite different: the 
south-eastern lobe is redder, to the point of becoming prominent in the near-IR, as revealed 



by the HST/NICMOS images of McCaughrean et al. (1998). It is not possible to explain 
such an inversion by assuming a standard flared disk model and common optical properties 
(scattering/absorption efficiencies) for the dust grains. However, a binary system with two 
stars having comparable luminosity and different effective temperatures could represent a 
more viable solution. Direct spectroscopy of the two lobes can address this possibility. 

Stellar encounters may also perturb circumstellar disks leading to the fomation of warps, 



spiral waves, tidal tails and disk truncation (Zinnecker & Yorke 2007). Those disturbances 



observed e.g. in RW Aur by Cabrit et al. (2006), are expected to occur in the ONC. Several 



studies have pointed out that in the dense inner regions of young clusters like the ONC stellar 
enconters may be frequent enough to produce major effects on circumstellar disks (see e.g. j 



Olczak, Pfalzner & Eckart 2008 and reference therein). 



Warps induced in the outer disk regions by the fly-by of another star propagate inward 



to a radius r < with a wavelength given by (Lubov, Ogilvie & Pringle 2002) 



c s (r) 
2VD 



(7) 



where c s (r) is the local disk sound speed and V and D are the relative velocity and closest 
approach distance of the encounter, respectively. The warp propagate inward to a radius 
r crit at which r = A(r), where the wavelength of the warp becomes larger than the disk 
size and gets reflected. Angular momentum, however, is still communicated to the inner 
disk so the disk becomes tilted, but not warped. The decay time-scale for the warp is given 
approximately by 
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where a is the dimensionless viscosity parameter, typically in the range a ~ 10- 2 -10- 4 (|HartrJiann 



2009). For our disk with = 400AU, assuming M+ = 1 M , the warp may be expected 
to last for a time of around 1.3 6 (a/10~ 4 )yr. Therefore, for the low values of the viscosity 
parameter a typical of the outer disk regions, a warp imprinted to the disk by a past close 
encounter may still be propagating at these later times. 

An interplay between the warped outer disk, the flat and tilted inner disk, the energy 
input from the central star(s) and the nebular environment supporting a conspicuous pho- 
toevaporative flow may account for the complex phenomenology observed in this unique 
system. In the near future, ALMA observations should greatly help clarifying the nature of 
114-426 by mapping with high spatial resolution the spatial distribution of the disk and its 
kinematic status. 



6. Conclusions 

We have used a set of HST images taken with WFPC2 and ACS to describe the complex 
morphology of the most prominent dark silhouette disk in the Orion Nebula, 114-426. The 
disk appears tilted, asymmetric and warped, with clear signs of photoevaporation of neutral 
material. Mapping the distribution of dust grains at the northern translucent edge of the 
disk using the Mie theory for standard circumstellar disk grains, we have found evidence for 
a spatial gradient of the grain radius, from 0.2 to 0.7 /im, with grains becoming smaller as 
the distance from the disk increases. These values are slightly below the minimal 1.9 /im 
grain size estimated by Shuping et al. (2003), but still compatible considering the difficulty 
of observing this source at infrared wavelengths. Overall, the disk morphology and location 
in the nebula indicate that this system is photoevaporated by the diffuse radiation field of 
the Orion Nebula, while being shielded from the radiation coming directly from the central 
Trapezium stars. Using our derived grain properties, we have estimated the mass loss rate 
from the disk surface and the time-scale for total disk dissipation, which turns out to be of 
the order of 10 4 yr. This short time scale, of the order of 1/100 of the cluster age, seems 
compatible with the exceptional nature of this system, namely its huge size and low mass. 
Finally, we briefly discuss the viability of possible mechanisms that may lead to the peculiar 
morphology of this system: external UV flux, binary star and past close encounter. 
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